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Abstract 



We generalize the nonlinear one-dimensional equation for a fluid layer surface 
to any geometry and we introduce a new infinite order differential equation 
for its traveling solitary waves solutions. This equation can be written as 
a finite-difference expression, with a general solution that is a power series 
expansion with coefficients satisfying a nonlinear recursion relation. In the 
limit of long and shallow water, we recover the Korteweg-de Vries equation 
together with its single-soliton solution. 

I. INTRODUCTION 

Shape deformations are important for an understanding of diverse many-body systems 
like the dynamics of suspended liquid droplets [1], long- lying excitations of atomic nuclei 
or "rotation- vibration" excitations of deformed nuclei and their fission or cluster emission 
modes [2] hydrodynamics of vortex patches [3], evolution of atmospheric plasma clouds [4], 
formation of patterns in magnetic fluids and superconductors [5], electronic droplets and 
quantum Hall effect [6], as well as resonant formation of symmetric vortex waves [7], etc. 

The theoretical description of such systems is often realized in terms of collective modes, 
such as large amplitude collective oscillations in nuclei, sound waves in solids, collective exci- 
tations in BEC, plasmons in charged systems, or surface moves in tubifex worms populations. 
Collective modes are especially important when their energies are lower than competing 
singe-particle degrees of freedom. Sometimes, however, single-particle or collective modes 
in the bulk of a system show particularly dense or sparse spectra. Systems of latter type 
are often referred as incompressible [1,5-7]. The incompressibility can be hard or relatively 
soft and still serves as a convenient limit for describing large differences in relevant length 
or time scales, as for the macroscopic motion of a liquid drop. Under these conditions, one 
can usually focus attention on the motion of the boundaries of the system, which usually 
has softer modes with frequencies that are lower than those of the bulk (e.g. Rayleigh waves 
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for a solid or surface waves for a liquid drop travel at speeds considerable slower than the 
bulk sound waves). 

Concentrating on the motion of the boundary of the system has considerable advantage: 
a reduction in the dimensionality of the problem that leads to relatively simple analytical 
treatments or a tremendous reduction in the effort needed to numerically solve or simulate 
the problem. Furthermore, associated with the incompressibility one usually finds micro- 
scopic conservation laws that expand into global constrains on the whole systems, even when 
the microscopic dynamics is completely local, as for example when the volume of a liquid 
drop is conserved. 

These shape deformations and their dynamics have played an important role in gaining a 
better understanding of numerous problems in diverse fields of physics. The incompressibility 
is normally reflected in the existence of a field that is piecewise constant, so that there is 
a sharp boundary between two or more distinct regions of the space with different physical 
properties. This field can be of classical origin like the density of a liquid or the charge 
density of a plasma, or it can originate in the quantum-mechanical properties of the system, 
like the magnetization of a type-II superconductor. Among the various examples where these 
questions are relevant, the waves on the surface of a liquid represent a unique opportunity to 
study the dynamics of shape deformations in a clean and controlled environment. We showed 
recently [1] that the surface dynamics of an incompressible liquid drop can be modeled by 
the Korteweg-de Vries (KdV) or the modified KdV (mKdV) equations for the circular shape 
[8,9]. The KdV equation was proposed a century ago, [9], to explain the dynamics of water 
waves observed in a channel by Scott Russel in 1834. 

The main property of this equation is the near-equal importance of nonlinear and disper- 
sion effects. Some of its most important solutions, like solitary waves, are single, localized, 
and non-dispersive structures that have a localized finite energy density. Among these solu- 
tions, the solitons are the solitary waves with special added requirements concerning their 
behavior at infinity (x — > oo, t — > oo) and having special properties associated with scat- 
tering with other such single-solitary wave solutions. Despite this nonlinearity, the KdV 
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equation is an infinite dimensional Hamiltonian system [10] and, remarkably, when the KdV 
solution evolves in time, the eigenvalues of the associated Sturm-Liouville operator d 2 + r\ 
(where 77 is the solution of the KdV equation) remain constant. Several extended and 
complex methods have been developed to study and to solve the KdV equation and other 
nonlinear wave equations like the nonlinear Schrodinger equation (NLS), the sine-Gordon 
equation (SG), etc. Among these there is the inverse scattering theory (1ST), theoretical 
group methods, numerical approaches. For a review of some of these techniques one can 
look into the book by Witham [11] or the paper by Scott et al [12]. 

Due to its properties the KdV equation, and its quadratic extension mKdV, were the 
source of many applications and results in a large area of non-linear physics (for a recent 
review see [13] and the references herein). However, in all these applications, the KdV equa- 
tion and its various generalizations (mKdV, KdV hierarchy, KP equation, supersymmetric 
generalizations, etc.) appear as the consequence of certain simplifications of the physical 
systems, especially due to some perturbation techniques involved. 

In the present paper we use only one of the two well-known necessary conditions for 
obtaining the KdV equation from a one-dimensional shallow water channel, i.e. only the 
smallness of the amplitude of the soliton r] with respect to the depth of the channel, h. 
This is the single condition which we use, with h taken to be an arbitrary parameter. In 
order to obtain the KdV equation from shallow water channel models, a second condition 
is imposed, i.e. the depth h of the channel should be smaller than the half-width A of 
the solitary wave. In the following, in our generalized approach, we do not impose this 
second condition nor the assumption of an infinite length channel; the channel length L, 
is taken to be an arbitrary parameter. By doing this we are actually investigating exact 
solutions for nonlinear modes in finite domains. What we actually do is replace the fixed- 
length boundary condition by a periodic requirement and in this way model general surface 
dynamics equations for liquid drops or liquid drop-like systems. This generalizations is 
important in connection with atomic and nuclear physics applications where the height of 
the perturbation may be comparable to the atomic or nuclear radius (big clusters formation 
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on the surface of symmetric molecules or nuclear molecules) so the dynamics of a shallow fluid 
layer becomes inappropriate. Therefore, we study, as a first step, the non-linear dynamics 
of a fluid of arbitrary depth in a bounded domain. This different starting point leads us to 
a new type of equation which generalizes in some sense the KdV equation (higher order in 
the derivatives and higher order nonlinearity) and also reduces to it the shallow liquid case. 

As we shall show, this generalized result (infinite order differential equation and higher 
order nonlinearity) can also be written in terms of a finite-difference equation. In this context 
we introduce here another possible physical interpretation for the translation operator in 
fluid dynamics, by relating its associate parameter to the depth of the fluid layer. In the 
present model we take into account not only the existence of a uniform force field (like 
an electric field or the gravitational field) but also the influence of the surface pressure 
acting on the free surface of the fluid. This implies that the KdV-like structure for the 
dynamically equation arrises in first order in the smallness parameter and one does not have 
to rely on a second order effect to see its consequences on the dynamics. We also note that, 
despite the general tendency followed by papers concerning applications of differential-finite- 
difference equations in physics of first introducing the equation and then searching possible 
applications, we obtained an infinite order partial differential equation and its connections 
to a finite-difference equation, in a natural way, starting from a traditionally physical one- 
dimensional hydro dynamic model. 

II. THE GENERALIZED KDV EQUATION 

Consider a one-dimensional ideal incompressible fluid layer with depth h and constant 
density p, in an uniform force field. We suppose irrotational motion and consequently the 
velocity field is obtained from a potential function y, t), that is V(x, y, t) = V$(:r, y, t) 
where we denote the two components of the velocity field V = (u, v). The continuity equation 
for the fluid results in the Laplace equation for $, AQ(x,y,t) = 0. This Laplace equation 
should be solved with appropriate boundary conditions for the physical problem of interest. 
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We take a 2- dimensional domain: x G [x — L, x + L] (as the "horizontal" coordinate) and 
y G [0, £(x,t)] (as the "vertical" coordinate), where x is an arbitrary parameter, L is an 
arbitrary length, and t) is the shape of the free surface of the fluid. The boundary 
conditions on the lateral walls x = xo ± L and on the bottom y — consist of the vanishing 
of the normal component of the velocity. The free surface fulfills the kinematic condition 

v|e = (6 + 6«)|e, (1) 

where we denote by £ the free surface of equation y = £(x,t) and the subscript indicates 
the derivative. Eq. (1) expresses the fact that the fluid particles which belong to the surface 
remain on the surface during time evolution. By taking into account the above boundary 
conditions on the lateral walls and on the bottom, we can write the potential of the velocities 
in the form: 

$(x, y,t) = jr, cosh (a k (t) cos (^^j + fa{t) sin (^^J j , (2) 

where a k and fa are time dependent coefficients fulfilling the condition 

. kn(x ±L) kn{x Q ±L) 
a k sin = fa cos , 

for any positive integer k. This restriction introduces a special time dependence of a k and 
fa, i.e. ^ — Ik =constant for any k, (3 ^ 0. If fa — we simply equate with the inverse 
of the above fraction. The a k (t) and fa{t) functions also depend on x . This special time 
dependence doesn't affect the general nature of the potential $, but does affect the balance 
between the two terms in the RHS of eq. (2). If we set all fa = 0, we can take arbitrary 
values for xq and general form for a k . In the infinite channel limit, L — > oo, there are no 
more restrictions concerning a k and fa functions, and we can choose x = without loss of 
generality. We introduce the function: 

f(x,t) = £ (-a fc (t) sin ^tt|) + fa{t) cos (^|)) ^, (3) 

so the velocity field can be written like: 
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u = $ x = cos(yd)f(x,t) 

v = $ y = -sm(yd)f(x,t). (4) 

where, for simplicity, the operator d represents the partial derivative with respect to the x 
coordinate. Equations (4) do not depend on L and therefore any approach toward the long 
channel limit must include the L — > oo (unbounded) limit. In the following we describe 
excitations of small height compared to the depth, and not necessarily large widths. Also, 
in the boundary condition eq.(l)we use velocities evaluated at y — £(x, t) — h + rj(x,t) to 
the first order in rj. Eqs.(4) reads: 

u{x, £(x, t),t) = [cos(hd) — r](x, t)dsm(hd)] f(x, t) 

v(x, £(x, t),t) — — [sm(hd) + r](x, t)d cos(hd)} f(x, t). (5) 

The dynamics of the fluid is described by the Euler equation at the free surface. The equation 
that results is written on the surface £ it in terms of the potential and differentiated with 
respect to x. By imposing the condition y — £(x, t), and by using a constant force field we 
obtain the form 

u t + uu x + vv x + gr] x + -^P x = 0, (6) 

where g represents the the force field constant and P is the surface pressure. Following the 
same approach as used in the calculation of surface capillary waves [17], we have for our 
one-dimensional case 

P l s = 15 = 7, aVx ^3/2 ~ ~ a Vxx, for small r/, (7) 
K (1 + Vx) 

where TZ is the local radius of curvature of the surface (in this case, the curvature radius 
of the curve y = £(x,t)) and a is the surface pressure coefficient. Inside the fluid the 
pressure is given by the Euler equation. The nonlinearities appear in the dynamics through 
the nonlinear terms in eqs. (1), (5), (6) and (7). Consequently, we have a system of two 
differential equations (1 and 6) for the two unknown functions: f(x,t) and rj(x,t), with u 



and v depending on 77 and / from eqs.(5). With / and rj determined and introduced in the 
expression for u and v,we can finally find the coefficients a k and f3 k . In the following we treat 
these equations in the approximation of small perturbations of the surface £, with respect 
to the depth, a = max\r]( k \x,t)\ « h, where k — 0, ...,3 are the orders of differentiation. 
In the linear approximation eqs. (1) and (6) become, respectively 

- sin(hd)f = rj t 

cos(hd)f t = -grj x + -r) xxx , (8) 
P 

and we obtain, by eliminating £ from eqs. (8) 

cos(hd)r] tt = sin(hd) (j^Vx - %xxx^j , (9) 

which, in the lowest order of approximation in (hd) for the sine and cosine functions, and in 
the absence of the surface pressure, gives us the familiar wave equation r} tt = c^rj XXl where 
Co = \fgh is the sound wave velocity. By introducing the solution r\ = e 1 ^-"*) in the 
linearized eq.(9) we obtain a nonlinear dispersion relation 



c, 2 = #(l + ^)^P. (10) 
pgl kh 



In the limit of shallow waters we find for the dispersion relation, in the a 7^ case of an 
acoustical wave as well as the a >> pg surface capillary wave limit, J 1 = —k 4 . In the 
absence of the surface pressure (a = 0) the function / is given in this linear approximation, 
at least formally, by 

which in the limit of a shallow fluid has a particular solution of the form 

f(x,t) = C ^ V . (12) 
For the time derivative of / we have, from the second equation of eq.(8) 

f t m (x, t) = (cos (hd))- 1 (-^rj x + % xxx \ , (13) 



which in the limit of a shallow fluid reduces to 



ft i X 't) — ^ Tj x + ^f]xxx- 

Following [10] we look for the solution of eqs.(l and 6) in the form 



a _ /a N 2 



f= h C ^ + {h) /2 ' 

/* = ~ C l\ ( COS ( hd )) 1 f > x + ~ ( COS (^)) 1 Vxxx + (j^j 92, 

which represents a sort of perturbation technique in |, where rj = afj. Of course a functional 
connection exists between the perturbation f 2 (x,t) and g 2 (x,t). Eq.(l) yelds in the lowest 
order in ^ 

— c sin (hd) fj = hfj t + ac (fj x cos (hd) 77 + 77 cos (hd) fj x ). (14) 

If we approximate sm(hd) hd — \ (hd) 3 , cos(hd) ss 1 — | (/i9) 2 , we obtain from eq.(14) 
the polynomial differential equation 

a?7i + 2c e 2 hr]r] x + Coe/ir/,, - c e— 77^ - — (VxVxx+VVxxx) = 0, (15) 

where e = f • The first four terms in eq.(15) correspond to the zeroth approximation for terms 
in /,obtained from boundary conditions at the free surface, in eq.(6. 1.15a) from Lamb's book 
[10], i.e., the traditional way of obtaining the KdV equation in shallow channels. In this case 
all the terms are first and second order in e. If we apply the second restriction with respect to 
the solution, that is the half-width is much larger than h, we can neglect the last parenthesis 
in eq.(15) and we obtain exactly the KdV equation for the free surface boundary conditions. 
In other words we understand the condition hd is "small" like (hd) f (x,t) « 1 over the 
entire domain of definition off. This means that the spatial extension of the perturbation 
f(x, t) is large compared to h, which is exactly the case in which the KdV equation arrises 
from the shallow water model (see Chapter 6 in [3], hdf(x, t) of order ps ^ = 5 « 1). 

By using again the approximations given by eq. (5) we can write the Euler equation eq. (6) 
in the form 
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mf + tydnf) + nf(d t Qf) + wf{mf) = -gvx + (16) 

where we use the notation 

Q = cos(hd) -r}dsin(hd), (17) 
oj = sin(hd) + r]d cos(hd) . 

Note that the operators given in eq.(17) satisfy the following interesting relation 

n + itu = e td \x + 2 (hd) + 2 (hri). (18) 

Following the same procedure as for the free surface boundary condition eq. (3) namely using 
an approximation for small i] we obtain from eq.(16) 

cos(hd)f t = -grj x + -rj xxx , (19) 
P 

which, in the lowest order in | and by using eq.(12 and 13) reduces to an identity. 

Before further analysis we would like to note that in the shallow water case, following the 
same notation as in Chapter 6 of Lamb [10], that is e = |, 5 = j, where / gives the order 
of magnitude of the half-width of the perturbation rj, and introducing a new parameter 
a = —-^r-p-, w e obtain for the Euler equation the form 

3 c 

fit' + Vx> + ~VVx' + u-fix'x'x' = 0, (20) 

which is again the KdV equation. The primes attached to the subscripts denote dimen- 
sionless units [10]. The difference between eq.(20) and the corresponding eq.(1.15.b) from 
Lamb's book is given by the inclusion of the surface pressure effects. If the coefficient a 
is large enough one can ignore second order terms, like for example 5 2 , to obtain the KdV 
equation. Of course, the above approach introduces changes in the differential equations 
which involve higher order perturbations like and from [10] or f 2 , g 2 in our case. 
The reduction of eq.(l) to KdV equation occurs if in eq.(14) we limit the expression to terms 
that are at most third order in hd 
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h 2 2c 

Vt + CoVx ~ c —T] xxx + —VVx = 0. (21) 
o h 

In the following we shall investigate this generalized KdV equation (gKdV) obtained 
from eq.(14) by keeping all terms in sin and cos, namely 

r]t + sin (/id) 77 + -^(rj x cos(hd)rj + rj cos(hd)r] x ) = 0. (22) 

lb lb 

Eqs.(l) and (16) yield, to higher orders in | and in hd, corresponding differential equations 
for the functions f 2 and g 2 , but here we shall study only eq.(21). 



III. THE FINITE-DIFFERENCE FORM OF THE GKDV EQUATION 

In this section we limit ourselves to the steady-state translational waves and consider 
only solutions of the form r](x, t) = r](x + Ac t) = r](X) where A G 1Z and X = x + Ac t. 
Eq.(21) can be written in the form 

Ahr} X {X) + v(X + *)-mX-*) + + ih) +, { X - ih) 

+ r,(xf {X + ' h)+ 1 (X -' h) =0 , (23) 
2 

if we suppose that i] is an analytic function in the domain Re(z) G (—00,00), ha(z) G 
(—h, h), z = x ± ih of the complex plane. We study rapidly decreasing solutions at infinity 
and we make the substitution v = e Bx for x G (—00, 0) and v = e~ Bx for x G (0, 00), where 
B is a positive constant. By introducing rj(X) = —hA + f{v) we obtain a differential-finite- 
difference equation (DFDE) for the function f(v) 

where we define the finite-difference operator as 

f(e iBh v) - fje^v) 

= e iBh v _ e - t Bk y • ( 25 ) 

We can write the solution of eq.(23) (or eq.(24)) as a power series in v 
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f(v) = J2 a n v n , (26) 

n=0 

and we choose ao = hA in order to have lim^-too rf(x) = 0. Equation (25) results in a 
non- linear recursion relation for the coefficients a n , that is 



Ahk + ^ sin (Bhk) 
B 



k-l 



= ~J2 n ( COS ( Bh ( k ~ n )) + COS ( Bh ( k ~~ a nClk-n- (27) 

n=l 

By taking k — 1 in the above relation, we obtain a x A/i + ^ sin (-B/i)] = 0. Without loss of 
generality 

and because of the arbitraryness of B we can write 

sin (Bh) 

A = (28) 

This relation fixes the velocity of the envelope of the perturbation if its asymptotic behavior 
is fulfilled. In order to have A ^ 0, we need Bh ^ kix for k integer. In this condition a\ is 
still arbitrary and by writing a k = a.k&i we have ct\ = 1 and the recursion relation 

2Bcos^^ *± Bh(2k-n-l) 

oth = 7— — -t 77— > ncos a n a.h- n , (29) 

jfc sin (Bh) - sin (kBh) ^ 2 n ' v ; 

for > 2. This recursion relation gives the coefficient for k in terms of those for k — 1 and 

lesser values. For a smooth behavior of the solution rj(X) at X = 0, i.e. continuity of its 

derivative, we must introduce the condition 

oo 

f v (l) = J2na n a n 1 ~ 1 = 0, (30) 

71=1 

or require that the derivative of the power series f(v) with coefficients given in eq.(28) to 
be zero in z G R, z — a±. This set the value for a±. 

In the following we study a limiting case of the relation eq.(28), by replacing the sin and 
cos expressions with their lowest nonvanishing terms in their power expansions 

6 k ~ l 

ftfc = wfc(F-i) g; ft " afe - (31) 
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It is straightforward exercise to prove that 

a * = (M^ k (32) 

is a solution of the recursion equation. This can be done using mathematical induction and 
by taking into account the relations 

k ^ 2 _ k(k- l)(2jfe-l) 

k -\ 3 (k(k-i)V 



E 



n 



n=l \ " J 

We can write the power expansion 

oo , i s k— 1 

which has the radius of convergence 72. = 2B 2 h 3 (due to the Cauchy-Hadamard criteria). 
The function g(z) can be written in the form 

g(z) = z(—± J -) 2B 2 h 3 = — Z -— 2 . (34) 

Condition eq.(31) results in a\ = —2B 2 h 3 and 

which provides 

oo , 

r ] (x) = 2B 2 h 3 Y / k(-e- B \ x \) 
k=i 

e -B\X\ 

2B 2 h 3 



[l + e- B \ x \) 2 



B 2 h 3 



2 • 



(36) 



2 Hh(^; 

As expected, this solution is exactly the single-soliton solution of the KdV equation and it 
was indeed obtained by assuming h small in the recursion relation eq.(28). So, the gKdV 
equation has two general features: in the reduction given by eq.(16) it yields the KdV 
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equation, and for eq.(230, namely in the limit hd "small" the differential equation and one 
of its solutions, r](X) go into the KdV equation and its single-soliton solution. 

In general we do not obtain a simple solution like the relation (33) but we have all the 
necessary information from eq.(28). It seems that the power series g(z) with the coefficients 
given by the recursion relation (28) has a nonvanishing radius of convergence. The problem 
of the existence of a real point z Q in the disk of convergence or on its edge (g'(z ) = 0) needs 
further study. For the KdV equation, this point is on the edge of the disk of convergence 
of the power series with coefficients given by the recursion relation (28). We mention that 
g'(z ) = means that the function is not univalent (it is not injective in the neighborhood 
of this point). 

An interesting observation concerning a conjecture formulated by Bieberbach in 1916 
seems to be worth noting: if f(z) = z + a 2 z 2 + . . . is analytic and univalent in the unit 
disk, then \a n \ < n for all n, with equality occuring only for rotations of the Koebe function 
k(z) = -. rf s x [14]. The form of our solution for the KdV equation, eq.(34), is exactly the 

(1 z) 

Koebe function up to a scaling of the variable z. This conjecture was proved in 1986 by de 
Branges [15]. This fact suggests that one could obtain information about the possible zeros 
of the g'{z) if it is possible to obtain estimates for the radius of convergence of the power 
series with coefficients given by eq.(30). 

The change of the variable that we used to obtain the nonlinear DFDE eq.(25) can 
be used in the time-dependent equation obtained from eq.(23). If r](x,t) = f(v,t) where 
v = e Bx or v = e~ Bx , then we obtain 

±Br(v) Sf J"^^ ± 2 sin (Bh) 5F(v, t) = 0, (37) 

where ± refers to the positive or negative real semiaxis. If we try a solution of the form 
o cifv k , we obtain the following infinite system of ODE: 
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Bh(n-l)X x , , 3 , , f(2k-n-l)Bh\ 
T 2Scos — A- ^ cos ^ ^-L J . (38) 

We note that if ao(to) — this relation is preserved in time. The above system can be 
resolved step-by-step when the equation for a n involves only the coefficients au with k < n. 
However, the smooth behavior of the solution in x = at any moment of time is a nontrivial 
problem. 

IV. CONCLUSIONS 

We have shown that the KdV equation used to describe fluids in shallow channels can be 
generalize to liquids flowing at any depth or length. As a consequence the role of the depth 
parameter in various applications that use such a model can be explored, like for example 
the dynamics of cluster formation and fission. The present paper announces two key results: 
a generalization of the KdV equation starting from a physical model; and anembeding of 
this nonlinear equation into a DFDE. 

By using a nonlinear hydrodynamic approach for a fluid layer of arbitrary depth and 
length, we obtained a differential equation of infinite order which is the generalization of the 
KdV equation for the corresponding shallow water pattern. We succeeded in rewriting this 
equationin the form of an finite-difference multiscale formalism, and consequently we have 
obtained a nonlinear recursion relation for the coefficients of its general traveling solution. 
We stress the importance of the introduction of the surface pressure term, eqs.(8,9) which 
yields the dispersion term, that is the term proportional to r] xxx m the equation. This term 
is essential in two respects: first it introduces the dispersion in a lower order of magnitude 
than in the traditional case [1], and second, it is responsible for dispersion in the case of 
the cylindrical geometry [10]. As it should be, the generalized KdV equation and its formal 
general traveling solution approach the KdV ones, in the shallow layer limit of the model. 

Based on these results we conjectured that there may exists a deep connection between 
NPDE, infinite order linear ODE and finite-difference equation, together with their symme- 
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tries. In short, we have shown that by starting from a one-dimensional model for an ideal 
incompressible irrotational fluid layer, a more general differential structure thanthe KdV 
equation is obtained and it has interesting properties, not all explored, that reduce as they 
approach the KdV equation in the shallow layer limit. 

We also stress that the result of section 3, namely the realization of the gKdV equation in 
terms of a DFDE could be the starting point for searching for more interesting symmetries. 
Finite-difference equations are related to self-similarity and wavelet bases, too. The B 
parameter in our model is similar with the scale parameter in multiresolution analysis, and 
we can obtain in this way a connection with this field. 

We anticipate many possible applications of such a formalism, especially in the field of 
clusters, fusion, fission, drops, bubbles and shells. If the gKdV equation could be shown to 
arise from a Lagrangean, or Hamiltonian formalism, one can apply this result in physical 
models which involve nonlinear shapes or contours. 

REFERENCES 

[1] A. Ludu and J. P. Draayer, Phys. Rev. Lett. 80 (1998) 2125; R. G. Holt and E. 
H. Trinh, Phys. Rev. Lett. 77 (1996) 1274; R. E. Apfel et al, Phys. Rev. Lett. 78 (1997) 
1912. 

[2] A. Bohr, Fys. Medd. K. Dan. Vidensk. Selsk, 26 (1952) 14; R. Gherghescu, A. 
Ludu and J. P. Draayer, J. Phys. G: Nucl. Part. 27 (2001) 63. 

[3] N. J. Zabusky et al, J. Comput. Phys. 30 (1979) 96. 

[4] E. A. Overman and N. J. Zabusky, Phys. Rev. Lett. 45 (1980) 1693. 

[5] A. T. Dorsey and R. E. Goldstein, Phys. Rev. B 57 (1998) 3059 

[6] C. Wexler and A. T. Dorsey, Phys. Rev. Lett, 82 (1999) 620, Phys. Rev. B 60 
(1999) 10971. 

[7] L. Friedland and A. G. Shagalov, Phys. Rev. Lett, 85 (2000) 2941. 
[8] M. Remoissenet, Waves Called Solitons (Springer- Verlag, New York, 1999); Y. 
Tsukuhara, et al, Appl. Phys. Lett, 77 (2000) 2926. 

16 



[9] D. J. Kortweg and G. de Vries, Phil. Mag., 39 (1895) 422. 

[10] R. Abraham and J. E. Marsden, Foundations of Mechanics (The Be- 
jamin/Cummings Publishing Company, Inc., Reading, Massachusets, 1978); G. L. Lamb, 
Elements of Soliton Theory (John Wiley & Sons, New York, 1980). 

[11] G. B. Witham, Linear and Nonlinear Waves (John Wiley & Sons, New York, 1974). 

[12] A. C. Scott, F. Y. F. Chiu and W. Mclaughlin, Proc. IEEE 61 (1973) 1443. 

[13] S. V. Vladimirov, M. Y. Yu and V. N. Tsytovich, Phys. Rep. 241 (1994) 1. 

[14] L. Bieberbach, S. -B. Preuss. Acad. Wiss. (1916) 940; K. Loewner, Math.Ann., 89 
(1923) 103. 

[15] L. de Branges, Acta Math., 154 (1985) 137. 



17 



